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Abstract 



We briefly review various mappings of fermion pairs to bosons, including 
^ . those based on mapping operators, such as Belyaev-Zelevinskii, and those 



< 

in 



on mapping states, such as Marumori; in particular we consider the work of 
Otsuka-Arima-Iachello, aimed at deriving the Interacting Boson Model. We 
then give a rigorous and unified description of state-mapping procedures which 



. allows one to systematically go beyond Otsuka-Arima-Iachello and related 

00 , 

. approaches, along with several exact results. 
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I. INTRODUCTION 

Professor Belyaev (with V. G. Zelevinskii) (BZ) pioneered the mapping of fermion sys- 
^ ; tems onto bosons more than thirty years ago These original attempts at bosonization 
I of fermion systems was motivated by collective particle-hole modes in nuclei. Since that 
time the interacting boson model 0] (IBM) has been phenomenologically very successful 
in explaining low energy nuclear spectroscopy for heavy nuclei. The bosons in this model 
are thought to represent monopole (J=0), quadrupole (J=2), and sometimes hexadecapole 
(J=4) correlated pairs of valence nucleons in the shell model. The IBM Hamiltonian is Her- 
mitian, usually has at most two-boson interactions, and conserves boson number, reflecting 
the particle-particle, rather than particle-hole, nature of the underlying fermion pairs. While 
one can numerically diagonalize the general IBM Hamiltonian, one of the strengths of IBM 
are the algebraic limits corresponding to the subgroups SU(3), U(5), or 0(6), with analytic 
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expressions for excitation bands and transition strengths, which encompass an enormous 
amount of nuclear data. 

The microscopic reasons for the success of such a simple model are elusive. Otsuka, 
Arima, and lachello, along with Talmi, have used a mapping of the shell model Hamiltonian 
to the IBM Hamiltonian |^,^ based on the seniority model but these attempts have not 
done well for well-deformed nuclei. For this reason we have revisited boson mappings to see 
if we can understand the success of the IBM starting from the shell model. 

In the next section we sketch out various historic approaches to boson mappings We 
follow Marumori and Otsuka et al. (OAI) in our mapping procedure which maps 
fermion states into boson states and construct boson operators that reproduce fermion ma- 
trix elements. We give the boson representation of the Hamiltonian and review the result 
that, in the full space, it factorizes into a boson image, which is the same as the BZ Hamilto- 
nian, times a normalization operator which projects out the spurious states. However, since 
our goal is to understand the IBM which only deals with a few of the enormous degrees of 
freedom of the shell model, we go on to discuss boson images in truncated spaces. This, we 
shall see, gives rigorous insight into the OAI mapping and shows how to go systematically 
beyond it. 



). As the fermion 



II. A BRIEF HISTORY OF BOSON MAPPINGS 

The fundamental goal is to solve the many-fermion Schrodinger equation 

i^l^A) =^a|^a) (1) 

and find transition matrix elements between eigenstates, ^aa' = (^^x T 
Fock space may be so large as to make direct solution intractable, the idea of a boson mapping 
is to replace the fermion operators with boson operators, using only a minimal number of 
boson degrees of freedom, that approximate the spectrum and transition matrix elements of 
the original fermion problem. There are two approaches to boson mappings which we now 
review. 



The first approach, epitomized by Belyaev and Zelevinskii |l| , is to map fermion operators 
to boson operators so as to preserve the original algebra. Specifically, consider a space with 
2^2 single-fermion states; gnify fermion creation and annihilation operators. The set 

of all bilinear femion operators, form the Lie algebra of 80(4^2), as embodied 

by the commutation relations 



d^dj , dj^dj^ 



[aittj, akai] = 
SiiSjk + Sikojaj + 6jialai - {i j) 



aittj, a\ai 
d^ dj J dj^di 



djkttiai - {i ^ j) 
- 5jka\ai — 5iia\aj 



(2) 
(3) 
(4) 
(5) 



At this point it is convenient to introduce collective fermion pair operators 

4-;^E(A9^^.ala]. (6) 

We always choose the VL{2VL — 1) matrices Ap to be antisymmetric to preserve the underlying 
fermion statistics, thus eliminating the need later on to distinguish between 'ideal' and 
'physical' bosons. We also assume the following normalization and completeness relations 
for the matrices: 



tr Ajj — 5o/3; 



(7) 
(8) 



Generic one- and two-body fermion operators we represent by T = Yl^ijTijajaj, V = 
J2iiu (/^l V |z/) Aj^A,^, where Tij = {i\ T from such operators one can construct a fermion 



Hamiltonian H. Now one has the following commutation relations: 



Aa,Ap 



0; 



Aa, A^ =Sai3-2j2 (aJjAo) . . aUj; 



A^,T\ =2^tr (A^TAUA 



fi,f2] =5:(TiT2-T2Ti),,aJa, 



(9) 
(10) 

(11) 
(12) 



The method of Belyaev and Zelevinskii is to find boson images of the bifermion operators, 
iK)B = &i + E <'b%h, + <''%blb\b,K + ... (13) 

0/37 a/B'jSe 

iA,)s = (4)L (14) 

a/3 alBjS 

where b^, 6^ are boson creation and annihilation operators, [b^, fcj^] = 5a/3, with the coefficients 
X, y chosen so that the images (Aj^) b, (A^) b, (T) b have the same commutation relations as in 
Because the algebra is exactly matched, if one builds boson states in exact analogy 
to the fermion states then the full boson Fock space is not spanned and one does not have 
nonphysical or spurious states. For the state-mapping methods described below, especially 
when the entire boson Fock space is used, identifying and decoupling from spurious states 
is an important and problematic issue. 

On the other hand, the Belyaev- Zelevinskii expansion is in general infinite. In the full 
boson Fock space, that is, no truncation of the boson degrees of freedom, the image of one 
body operators is finite and given quite simply by {T)b = 2 Sa/jtr {AaTAp)bl^bf3. Since any 
fermion Hamiltonian can be written in terms of one-body operators, the boson image of a 
finite fermion Hamiltonian can be finite in the full space. The states that one must use then 
are built from the pairs given in ([13|) which will not just be products of bosons but will 
include exchange terms. For example, for two bosons and using (0), 

AlAUo) ^ (bibU x'^J%bl) \0) . (16) 



These exchange terms are due to the antisymmetry. We shall take care of such exchange 
effects by introducing a norm operator in the boson space. For truncated spaces, however, 
the expansion of the BZ Hamiltonian is infinite. 

The Dyson mapping |p is a variant of Belyaev-Zelevenskii, in which one makes the 
mapping 

Aa^ba] (17) 



4 ^ fct _ 2 ^ tr (A,A^A,At )6i6t 5^ (18) 

f ^Y.'^ti{K^TA\)h%. (19) 

(20) 

The operators are then clearly finite; on the other hand they are just as clearly non- 
Hermitian. From a computational viewpoint non-Hermiticity is only a minor barrier, but it 
is an obstacle to an understanding of the microscopic origin of Hermitian IBM Hamiltonians. 
Furthermore the Dyson operators mix spurious and physical spaces. 

Marshalek points out there exist mappings that are both finite and Hermitian, but 
these in general require projection operators to eliminate spurious states. We will regain 
this result later on in this paper. 

The second major approach, pioneered by Marumori 0, is to map fermion states and 
construct boson operators that preserve matrix elements. For fermion many-body (shell- 
model) basis states one often uses Slater determinants, antisymmetrized products of single- 
fermion wavefunctions which we can write using Fock creation operators: ' ' ' |0) for 
n fermions. For an even number of fermions one instead constructs states from products of 
fermion pairs, 

1^/^) = n 4. lo) ; (21) 

m=l 

if the number of fermion is fixed at n then m runs from 1 to = The original work 

of Marumori, however, focused on particle-hole excitations and so the number of pairs and 
consequently bosons is not fixed. These states are not trivially orthonormal. They must be 
orthogonalized. Exactly how this orthogonalization is accomplished will be a key theme in 
this paper. 

Marumori constructs the norm matrix 

M^p = {^a\'^p) (22) 

and then the Usui operator 
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u=Y: \^,)i^f)-,'f{^a\ (23) 

where bosons states are constructed in strict analogy to the fermion states, 

N 

m= n^ftJO). (24) 

m=l 

Then the Marumori expansion of any fermion operator is 

Ob = UOpUl (25) 

Clearly Marumori is best suited for particle-hole states with only a few excitations. If one 
applied it to a system with numerous particle-particle pairs, as for the IBM, one obtains 
clumsy many-body terms. Kishimoto and Tamura addressed this last issued by intro- 
ducing a "linked-cluster" expansion which they then grafted into a BZ-type scheme. 

Otsuka, Arima, and lachello (OAI), along with Talmi 0,^, investigated the microscopic 
origins of the Interacting Boson Model through boson mappings. Although they also mapped 
states, they differed from Marumori in some key details. First of all, they built states built 
on a fixed number of particle-particle, not particle-hole, pairs, and restricted the pairs to 
one monopole (J'^ = 0+) and quadrupole (J'^ = 2+) pair. These states were orthogonalized 
based on seniority. That is, they construct, for 2N fermions, low-seniority basis states of S 
and D fermion pairs, S'^~"'*D"''^, and then orthonormalize the states such that the zero- 
seniority state is mapped to itself, and states of higher seniority v are orthogonalized against 
states of lower seniority, 

\v) I V) = \v) + \v-2) + \v-4:) + ... (26) 



Then OAI calculate the matrix elements ( \Hf\ ) for n^, n^^ = 0, 1, 2 

and obtain the coefficients for their one plus two-boson Hamiltonian. These coefficients 
have an explicit A^-dependence (and for large A^ and arbitrary systems such matrix elements 
are not trivial to calculate, especially in analytic form!) and thus implicitly a many-body 
dependence. At first sight this is not entirely unreasonable as it is well known the IBM 
parameters change substantially as a function of the number of bosons, even within a major 
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'shell'. Nonetheless the OAI mapping has three drawbacks. The first is that it's not clear 
how to systematically calculate many-body contributions beyond that contained in the OAI 
prescription, whereas the method we shall describe is fully and rigorously systematic. The 
second is that the OAI prescription can induce many-body effects where none are needed. 
This point will be illustrated in section |1V C\ Thirdly, only the = 0, 1, 2 space is exactly 
mapped, but very deformed systems will involve large Ud- In fact, for an axial rotor limit, 
the average number of d-bosons in the ground state band is 2/3 the total number of bosons. 
As an alternative to OAI, Skouras, van Isacker, and Nagarajan proposed a "demo- 



cratic" mapping where the orthogonalization is based on eigenvectors of the norm matrix 
rather than seniority. 

In what follows we attempt to rigorously unify all the state-mapping methods. We have 
three strong results. First, we give general expressions for fermion matrix elements via 
boson representations. Second, we show how in several cases one can have exact, finite, and 
Hermitian boson images of fermion operators. Finally, we show how to extend both the 
OAI and democratic mappings in a systematic and rigorous fashion, and illustrate how the 
choice of orthogonalization can affect the many-body dependence of the boson images. 

III. BOSON REPRESENTATIONS OF FERMION MATRIX ELEMENTS 

The starting point of any state-mapping method is the calculation of matrix elements of 



fermion operators between states constructed from fermion pairs of the form (|2T|) , including 
the overlap: (\I^a|^/3), (^a H (^'^a T \l//3^, and so on. These matrix elements are 

much more difficult to compute than the corresponding matrix elements between Slater 
determinants. As we shall show, however, full and careful attention paid to the problem of 
calculation matrix elements can yield powerful results. Silvestre-Brac and Piepenbring |]12 



laboriously using commutation relations, derived a Wick theorem for fermion pairs. Rowe, 
Song and Chen [|13| using 'vector coherent states' (we would say fermion-pair coherent states) 

/ -s \ Af 

found matrix elements between pair-condensate wavefunctions, states of the form [An |0) . 
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Using a theorem by Lang et al. [|1^], we have generahzed [|T5| the method of Rowe, Song and 
Chen and recovered (actually discovered independently) the expressions of Silvestre-Brac 
and Piepenbring. Specifically, we construct generating functionals by taking the matrix 
element 







exp 



,(e...),(e...i 







oo / 9^fc— 1 

E^^tr 

k=l 



k 



(27) 



By taking derivatives of e^, etc., one computes the desired matrix elements in analytic form 
||15||. For pair condensate wave functions one can calculate the matrix elements iteratively 



and propose a variational principle |]T3[. Such a variational principle would be useful in 
determining the "best" microscopic structure for a truncated set of pairs. Alternately, 
Otsuka and Yoshinaga |]l6l choose their S and D pairs from Hartree-Fock-Bogoliubov states; 
the two approaches can probably be related in some approximation. This may be important 
is answering a basic question of IBM, the origin of algebraic limits: do they arise from 
changes in pair structure, or from effective many-body effects, or both? 

We now want to translate the fermion matrix elements into boson space. We take the 
simple mapping of fermion states into boson states 



N 



(28) 



m=l 



where the are boson creation operators. We construct boson operators that preserve 
matrix elements, introducing boson operators Tb, Vb, and most importantly the norm 
operator Mb such that ($q|7b|$/3) = (^a|T|^/3), (^alVsl*/?) = (^q|V"|*/3). and 
($q,|A/b|$/3) = {^a l^/s)- We term Tb,Vb the boson representations of the fermion op- 



erators T, V. One finds the 'linked-cluster' (a la Kishimoto and Tamura |10| although with 



differences) expansion of the representations to be of the form |jT5 

oo e e 

e=2 {a,T} i=l j = i 



(29) 
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and similarly for Vb, T'b- In the norm operator the £-body terms embody the fact that the 
fermion-pair operators do not have exactly bosonic commutation relations, and act to enforce 
the Pauli principle. The coefficients etc. can be written in closed, albeit complicated, 
form [p!5| . 

The norm operator can be conveniently and compactly expressed [lT7| , p!5| , p!8| in terms 
of the kih order Casimir operators of the unitary group SU(2i7), Ck = S'^tr (P)'^, P = 
J2aT bj^brA^Al (and so is both a matrix and a boson operator; the trace is over the matrix 
indices and not the boson Fock space) 



/ 1 (-1)'' . 



(30) 



where the colons ':' refer to normal-ordering of the boson operators. This norm operator 
takes into account the exchange terms in the BZ expansion of a fermion pair given in (|13|). 
Similarly — and this is a new result we have not seen elsewhere in the literature — the 



representations Tb, Vb can also be written in compact form |ll5| , |l8 



= 2 tr [A^TAt g] blKJ^B: 
Vi? = E 1^1 ^) E- [A<xAt g] tr [A.At G 
+ 4tr [A^AiPGA^AtG]}6t&.Mj:, 



(31) 



(32) 



where G = (1 + 2P)^^. These compact forms are useful for formal manipulation. Further- 
more they have the powerful property of exactly expressing the fermion matrix elements 
under any truncation, a fact not previously appreciated in the literature even for the norm 
operator By this we mean the following: suppose we truncate our fermion Fock space 
to states constructed from a restricted set of pairs {a}. Such a truncation need not corre- 
spond to any subalgebra. Then the representations in the corresponding truncated boson 
space, which still exactly reproduce the fermion matrix elements and which we denote by 
[J^b]t etc., are the same as those given above, retaining only the 'allowed' bosons with 
unrenormalized coefficients. For example 
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[A^^]^=:exp(--;^^^[c4) : (33) 



2 

^ k=2 



where 

[C,]t = 2:tr ([P]^)'=: , [P]^ = ^^fe.A^At. (34) 

CfT 

This invariance of the coefficients under truncation will not hold true for the boson images 
introduced below. 

With the boson representations of fermion operators in hand, one can express the fermion 
Schrodinger equation ([^) with = T + as a generahzed boson eigenvalue equation, 

nB\<^x)=ExArB\<^x)- (35) 

Here 'Hb is the boson representation of the fermion Hamiltonian. Every physical fermion 
eigenstate in (|1|) has a corresponding eigenstate, with the same eigenvalue, in (pS]). Because 
the space of states constructed from pairs of fermions is overcomplete, there also exist 
spurious boson states that do not correspond to unique physical fermion states. These 
spurious states will have zero eigenvalues and so can be identified. The overcompleteness 
also means that (|35D is harder to solve exactly than (|I]). So one truncates the model space. 

IV. BOSON IMAGES 



In general the boson representations given in (p9D, ( |3TD and ( |32D do not have good 
convergence properties, so that simple termination of the series such as (E^ in i-hodj terms 



is impossible and use of the generalized eigenvalue equation (|35|) , as written, is problematic. 



Instead we "divide out" the norm operator to obtain the boson image, i.e. schematically, 

h ~ ''Hb/J^b" (36) 



That this is reasonable is suggested by the explicit forms of ( |3lD and (0). The hope of 
course is that h is finite or nearly so, so that a 1+2-body fermion Hamiltonian is mapped 
to an image 
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h ~ Oib^b + e2b^b^bb + Osb^b^b^bbb + O^b^b^b^b^bbbb + ... (37) 

with the £-body terms, £ > 2, zero or greatly suppressed. We now discuss how this "dividing 
out" is to be carried out. 



A. Exact results: Full Space 

It turns out that for a number of cases the image of the Hamiltonian is exactly finite. 
In particular, for the full boson Fock space the representations factor in a simple way: 
Tb = MbTb = TbMb and Vb = J^bVb = VbN'b-, where the factored operators Tb, Vb, 
which we term the boson images of T, V , have simple form 



fB = 2j2tT (A^TAt)b%, (38) 



(39) 



blK + 2 E E (A^At A., At A.aU blbl^br' 



This image Hamiltonian Hb = Tb + Vb is the one determined by BZ if one decomposes the 
Hamiltonian into multipole-multipole form and then maps these mulitpole operators. As 
discussed earlier, these BZ multipole operators are finite in the full space. This result, and 
its relation to other mappings, was noted by Marshalek 

Thus any boson representation of a Hamiltonian factorizes: Tis = MbHb in the full 
space. Since the norm operator is a function of the SU(2f2) Casimir operators it commutes 



with the boson images of fermion operators |lgJT^, and one can simultaneously diagonalize 
both TIb and Nb- Then Eqn. (|35D becomes 

Hb\^x) = K\^x). (40) 

where E'^^ = Ex for the physical states, but E'^ for the spurious states is no longer necessarily 
zero. The boson Hamiltonian Hb is by construction Hermitian and, if one starts with at 
most only two-body interactions between fermions, has at most two-body boson interactions. 
All physical eigenstates of the original fermion Hamiltonian will have counterparts in (^0|). 
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It should be clear that transition amplitudes between physical eigenstates will be preserved. 
Spurious states will also exist but, since the norm operator A/^ commutes with the boson 
image Hamiltonian Hb-, the physical eigenstates and the spurious states will not admix. 
Also the spurious states can be identified because, while they will no longer have zero 
energy eigenvalues, they will have eigenvalue zero with respect to the norm operator. 



B. Exact Results: Truncated space 



The boson Schrodinger equation (p]|), though finite, is not much use as the boson Fock 
space is still much larger than the original fermion Fock space, and we still must truncate 
the boson Fock space. Although the representations remain exact under truncation, the 



H 



B 



T 



B 



H 



B 



T 



This was 



factorization into the image does not persist in general: 
recognized by Marshalek (An alternate formulation 0] does not require the complete 
Fock space, but mixes physical and spurious states and so always requires a projection 
operator.) 

If the truncation scheme represents a closed subalgebra (specifically, if the truncated set 



of fermion pairs are closed under double commutations) then a factorization [11^ 



n 



B 



B 



h 



D 



(41) 



does exist, with h£, at most two-body, but not necessarily Hermitian. We term it a Dyson 
image @J^. Under more restricted conditions on the structure of the pairs and the Hamil- 
tonian one can guarantee %£> is Hermitian and commutes with As ^. In the full space, of 
course, all definitions of boson images coincide and yield the same result. 

First, consider a partition of the single fermion states labeled by i = {ia,ic), where the 
dimension of each subspace is 2Qa, so that Q = 2Qa^c- We denote the amplitudes for the 
truncated space as and assume they can be factored, (Ajj)jj = (K.^)i^j^ (g) (Ag)j^j^, with 

[—lyK, where p = (symmetric) or p = 1 (antisymmetric). 
Furthermore we assume the completeness relation M), which was crucial for proving that 



KtK = KK^ = 2^ and 



Hb = J^BtlB lll5| , [T8[| , is valid for the truncated space; i.e.. 
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E(aD..>(a.) 



1 



5j jf ov 



(42) 



The norm operator in the truncated space then becomes 

-tr(P^):, 



k=2 ^"a- 



k 



(43) 



where P = J2af &5^&f A^-A^ so that [P]^ = (^^^ P- In this case the boson image of a 
one-body operator is the truncation of the boson image in the full space, 



Tn 



(44) 



2^tr (A^TAt) blW. 



(45) 



The representation of a two-body interaction can be factored into a boson image times the 
truncated norm, 



Nb 



however, while finite (l-|-2-body), is not simply related to 
one-body operators. If one writes 



(46) 



as is the case for 



(47) 



cr,T 



then matrix elements of the two-boson interaction are 



{(7(7'\v\tt') — 

Q„(2Q„ - (-i)p)(n, + 

+2Q„[tre(A^At A^,At A,A{,) - tr,(A^At A^, At A.KA^,)] 

-^la{2^la + (-l)P)tre(A,KAt A^,A|,)5<., J. 

Upon inspection one sees the image is not constrained to be Hermitian. 

Consider the additional condition between the matrix elements of the interaction: 



(48) 



13 



E(HV'k)E(A.)..,,,.,(Ai)^^^.,._,, 

H'l'^ In . 1 11 I, ^ . I ^ 



(49) 

where the factor Na = Qai'^^a + (^1)^) is the number of pairs in the excluded subspace. 
While condition (]49| ) looks complicated there are interactions that satisfy it; for example, 
two-body interactions constructed from one-body operators V = T^fl"^,^, where T^^ = 



Ai,A0 . When 
simply related: 



is satisfied then vo is Hermitian and although vd 



V, 



B 



they are 



vd 



(T,r 



V, 



B 



(50) 

by a factor which ranges 



with Jq^ = AQ'^/Na renormalizing the two-boson part of 
from unity (full space) to 2 for a very small subspace. Not all interactions satisfy (^9|); 
for example, the pairing interaction never does except in the full space. For the pairing 

ypairing 



mteraction (ft 



V ) = ^tifi^iyfiG, and AqAJ = and the image (^6]) {jg^""™^ becomes 
(remembering Q = 2Qa^c) 

g|no[1-§N + ^ + ^]+ E tr (A^At AoA^,)4&SM.' 



n n 



(51) 



where N is the total number of bosons, N = J^fbfbf, and No = ^0^0 • The second term in 



the above is not Hermitian but can be transformed away by a similarity transformation [20 



leaving the first term as a finite Hermitian image which gives the correct eigenvalues for all 
N. 



The SO (8) and Sp(6) models |^ belong to a class of models which have a subspace 
for which (|4^ ) is valid and interactions which satisfy (|49|). In these models the shell model 
orbitals have a definite angular momentum j and are partitioned into a pseudo orbital angular 
momentum k and pseudospin i, j = k + i. The amplitudes are then given as products of 
Clebsch- Gordon coefficients, 
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i^V).. = 7y {kmi,kmj\K^M^) (52) 

where K and / are the total pseudo orbital angular momentum and pseudospin respectively 
of the pair of nucleons. For the SO (8) model i = | and one considers the subspace of pairs 
with K = {p = 0), {Al,)ij = ^^^-^^^'^^^{i fii,i fij\Ia fia)] in the Sp(6) model k = 1 and one 
considers the subspace with / = (p = 1), (^4^)^^ = (^+(~^) ^ (^kmi,kmj\KaMa). The 
complicated conditions p9|) hold true for important cases, for example, the quadrupole- 
quadrupole and other multipole-multipole interactions in the S0(8) and Sp(6) models (that 
is, interactions of the generic form ■ in the notation of ^T|) have Hermitian Dyson 
images. Not all interactions in these models have Hermitian Dyson images. For example, 
pairing in any model (see ( |^ )and, in the S0(8) model, the particular combination gQ{S'^S+ 
ip2 . p2-j -^j^icj^ ig ^j^g SO (7) limit. It so happens that these particular cases nonetheless can 
be brought into finite, Hermitian form as discussed in the next section. 



C. Approximate or numerical images 



The most general image Hamiltonian one can define is 



h = U 



B 



-1/2 



T 



n 



B 



Ms 



B 



-1/2 



(53) 



which is manifestly Hermitian for any truncation scheme and any interaction, with U a 
unitary operator. (Because the norm is a singular operator it cannot be inverted. Instead 



-1/2 . 



is calculated from the norm only in the physical subspace, with the zero eigenvalues 



which annihilate the spurious states retained. Then h does not mix physical and spurious 
states.) If W = 1, this is the democratic mapping Again, for the full space h = Jid = Hb- 
This prescription is, we argue, useful for a practical derivation of boson image Hamil- 
tonians. Ignoring for the moment the unitary transformation W, consider the expansion 



of h. The operators 



B 



and 



B 



-1/2 



have similar expansions, and by multiplying 



out (|53D one sees immediately that the coefficient 6i depends only on up to £-body terms 

-1/2 

, derived from 2£-fermion matrix elements which are tractable for £ 
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in [Hb\t and 



small. Ideally h would have at most two-body terms, and our success in finding finite images 
in the previous section gives us hope that the high-order many-body terms may be small; at 
any rate the convergence can be calculated and checked term-by-term. Specifically, consider 
the convergence of the series (^) as a function of £. A rough estimate is that, for an A^- 
boson Fock space, one can truncate to the £-body terms if for ^' > i, 6gi is sufficiently small 
compared to x (A^ — i')\/[N — the strictest condition is to require O^i <^ 0^/{l' — €)\. 



The Hermitian image /i, defined in (|53|), is related to the Dyson image /id, defined in 
(PI) , by a similarity transformation S =U Mb ^ , 

h = ShDS-\ (54) 

The similarity transformation S orthogonalizes the fermion states \^ a) inasmuch 
{S-^yUsS-'^ = 1 in the physical space (and = in the spurious space). This is akin 
to Gram-Schmidt orthogonalization and the freedom to choose U, and S, corresponds to the 
freedom one has in ordering the vectors to be Gram-Schmidt orthogonalized. The OAI and 
democratic mappings are just two particular choices out of many; the latter takes U = 1. We 
can use the freedom in the choice ofU to our advantage. Consider the S0(8) model PT| and 



its three algebraic limits: the pure pairing interaction, the quadrupole ■ interaction, 
which can be written in terms of SO (6) Casimir operators, and the linear combination of 
pairing and quadrupole S'^S + ■ P^ which can be written in terms of S0(7) Casimirs 
for details and notation). As discussed in the last Section, the Dyson image of the 



see 



quadrupole interaction is Hermitian and finite, and ho = h with U = 1. The Dyson images 
of the pairing and SO (7) interactions are finite but non- Hermitian. We have found Ws ^ 1 
for both these cases (but not the same U) such that their respective Hermitian images h are 
finite; the one for pairing is exactly the OAI prescription, while that for SO (7) is exactly 
opposite, orthogonalizing states of low seniority against states of higher seniority. These 
general Hermitian images do not have a simple relation to the truncated full image, as do 
the Hermitian Dyson image (|5D|). The one-body piece remains unchanged but there can be 
significant renormalization, and even change of sign, of the two-body piece. For example, for 
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the pairing interaction Hb ^ — s^s + -^s^s^ss-\- additional terms, including off-diagonal 
terms such as d^d'^ss, whereas h = s'^s — ^s^s'^ss+ (depending on the truncation) terms such 
as d^d^dd but no off-diagonal terms. Hence we see here the renormalization is not just a 
simple overall factor / as it was for the Hermitian Dyson image: it is — 2Q for the s^s^ss 
term but for terms such as d^d'^'ss. 
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FIGURES 



FIG. 1. Spectrum of SO (7) interaction, for 7 bosons, in SO (8) model with exact (left) and 
approximate (right) two-body boson Hamiltonians. 

If one uses an "inappropriate" transform S it can induce an unneeded and unwanted 
many-body dependence. This principle we illustrate in the SO (8) model with the SO (7) 
interaction, whose spectrum is exactly known and for which we can derive a finite Hermitian 
image with no many-body dependence; this is the left-hand spectrum in Figure 1. For 
the right-hand side we took the U appropriate for pairing, that is the OAI prescription 
determined for TV = 2, and calculated the spectrum for = 7 keeping only the strict two- 
body terms. The distortion in the spectrum from the exact result, such as the overall energy 
shift and the large perturbation in the third band, indicates missing many-body terms. That 
is, if one mapped the SO (7) interaction using the canonical OAI procedure one would find 
of necessity a many-body dependence in the interaction coefficients. By orthogonalizing the 
basis in the a different way, however, as expressed by a different choice of U, the many- 
body dependent vanishes. Therefore it is possible that some of the A^-dependence of OAI is 
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induced by their choice of orthogonahzation and could be minimized with a different choice. 
We are currently exploring how to exploit this freedom to best effect. 



V. SUMMARY 

In order to investigate rigorous foundations for the phenomenological Interacting Boson 
Model, we have presented a rigorous microscopic mapping of fermion pairs to bosons, paying 
special attention to exact mapping of matrix elements, Hermiticity, truncation of the model 
space, and many-body terms. First we presented new, general and compact forms for boson 
representations that preserve fermion matrix elements. We then considered the boson image 
Hamiltonian which results from "dividing out" the norm from the representation; in the full 
boson Fock space the image is always finite and Hermitian; in addition we discussed several 
analytic cases for truncated spaces where the image is also finite and Hermitian. Finally, 
we give a prescription which is a generalization of both the OAI and democratic mappings; 
in the most general case for truncated spaces the Hermitian image Hamiltonian may not 
be finite but we have demonstrated there is some freedom in the mapping that one could 
possibly exploit to minimize the many-body terms. This freedom, which manifests itself in a 
similarity transformation that orders the orthogonahzation of the underlying fermion basis, 
depends on the Hamiltonian. 

This research was supported by the U.S. Department of Energy. The boson calculations 



for Figure 1 were performed using the PHINT package of Scholten [22 
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